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Abstract. Remnants of neutron-star mergers are essentially massive, hot, differentially rotating neutron 
stars, which are initially strongly oscillating. As such they represent a unique probe for high-density matter 
because the oscillations are detectable via gravitational-wave measurements and are strongly dependent 
on the equation of state. The impact of the equation of state for instance is apparent in the frequency 
of the dominant oscillation mode of the remnant. For a fixed total binary mass a tight relation between 
the dominant postmerger oscillation frequency and the radii of nonrotating neutron stars exists. Inferring 
observationally the dominant postmerger frequency thus determines neutron star radii with high accuracy 
of the order of a few hundred meters. By considering symmetric and asymmetric binaries of the same 
chirp mass, we show that the knowledge of the binary mass ratio is not critical for this kind of radius 
measurements. We perform simulations which show that initial intrinsic neutron star rotation is unlikely 
to affect this method of constraining the high-density equation of state. We also summarize different pos¬ 
sibilities about how the postmerger gravitational-wave emission can be employed to deduce the maximum 
mass of nonrotating neutron stars. We clarify the nature of the three most prominent features of the post¬ 
merger gravitational-wave spectrum and argue that the merger remnant can be considered to be a single, 
isolated, self-gravitating object that can be described by concepts of asteroseismology. We sketch how 
the consideration of the strength of secondary gravitational-wave peaks leads to a classification scheme of 
the gravitational-wave emission and postmerger dynamics. The understanding of the different mechanisms 
shaping the gravitational-wave signal yields a physically motivated analytic model of the gravitational- 
wave emission, which may form the basis for template-based gravitational-wave data analysis. We explore 
the observational consequences of a scenario of two families of compact stars including hadronic and quark 
matter stars. We find that this scenario leaves a distinctive imprint on the postmerger gravitational-wave 
signal. In particular, a strong discontinuity in the dominant postmerger frequency as function of the total 
mass will be a strong indication for two families of compact stars. 

PACS. 26.60.Kp Equations of state of neutron-star matter - 97.60. Jd Neutron stars - 04.30.Db Wave gen¬ 
eration and sources - 95.85.Sz Gravitational radiation, magnetic fields, and other observations - 04.25.dk 
Numerical studies of other relativistic binaries - 21.65.Qr Quark matter 


1 Introduction 


becomes stronger, and thus the decay of the orbit pro¬ 
ceeds increasingly faster. For PSR B1913-I-16 the orbital 
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By means ot intense experimental ettorts gravitational waves j • o ms rm i- 

/A-tur.N ,. 1 , 1 . 1 . . .. . 1 -.1 ? .1 . . milliseconds in the next 3 x 10 years, ihe corresponding 

GW emission reflects the dynamics of the orbital motion 

and thus increases in frequency and amplitude resulting 


(GWs) are expected to be detected within the next years. 

To date, their existence is only indirectly proven by the ob¬ 
servations of neutron star (NS) binaries II191III41II27I . in . , . ... . , , • . • -iixiv 

tt 1 r,n 1 V- 1 m ni o : a chirp-like signal, which is determined by the binary 

particular the Hulse-Taylor binary pulsar PSR B1913-I-16 lDJ. . “ ’ ^ .... . . ™ -i 
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The orbital period and orbital separation of this binary 
decrease precisely as determined by General Relativity, 
which predicts that the GW emission continuously re¬ 
duces the system’s angular momentum and energy. As the 
orbital frequency increases with time the GW emission 
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and 6 in [H]). Only during the last seconds of this so- 
called inspiral (named after the shape of the stellar tra¬ 
jectories) the signal will enter the sensitivity band of the 
existing and upcoming GW detectors between roughly 10 
and 10000 Hz mmm- Importantly, the total popula¬ 
tion of NS binaries within a distance of a few 100 Mpc 
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to our Galaxy (i.e. the range that will be covered by GW 
instruments) is expected to be sufficiently high such that 
per year of the order of 40 binaries in the last phase of 
their inspiral will enter the sensitivity band of GW detec¬ 
tors [5]. 

The continuously decreasing orbital separation of NS 
binaries inevitably leads to the merging of the binary com¬ 
ponents. On a dynamical time scale of the order of a mil¬ 
lisecond the two stars form a single, massive, differentially 
rotating object (see e.g. Fig. Eland Fig. 2 in [12])- Mostly 
because of its rotation the merger remnant can be sup¬ 
ported against the gravitational collapse even if the total 
binary mass significantly exceeds the maximum mass of 
non-rotating NSs. Only for very high total binary masses 
a prompt gravitational collapse occurs. The stellar object 
forming during a coalescence of two NSs is very interest¬ 
ing from the perspective of high-density matter physictQ 
because it resembles a massive, hot, rotating NS, which is 
strongly oscillating (see e.g. [^[8l [T^[52ll 106) for reviews). 
The oscillations are induced by the merging process, no¬ 
tably the fundamental quadrupolar fluid mode of the rem¬ 
nant is strongly excited and dominates the postmerger 
GW signal. The frequencies of the different excited oscil¬ 
lation modes are characteristic of the remnant’s mass and 
of the high-density equation of state (EoS). The mass of 
the remnant is approximately identical to the total binary 
mass. The binary mass as well as (in principle) most of 
the oscillation frequencies of the merger remnant are ob- 
servationally accessible by GWs. Therefore, NS mergers 
offer the possibility of NS asteroseismology and of inferring 
properties of the high-density EoS, or equivalently, of stel¬ 
lar properties of non-rotating NSs, via observations with 
the upcoming and existing GW detectors. The merger rate 
is estimated to be roughly of the order of 10“® events per 
year per Milky way equivalent galaxy, which corresponds 
to a detection rate of about 40 inspiral detections per year 
for Advanced LIGO and Advanced Virgo when the instru¬ 
ments operate at their design sensitivity (see e.g. [2] for a 
collection of rate estimates and their uncertainties). With 
these instruments the oscillations during the postmerger 
phase are detectable for relatively nearby events [35l[36], 
which is an exciting prospect given that already a sin¬ 
gle event is sufficient to provide tight constraints on the 
EoS [niiii]- 

Because of its high mass, the postmerger remnant is 
particularly interesting since the central densities are rel¬ 
atively high compared to those of many observed NSs in 
binary systems, which have masses in the range of 1.2 Mq 
to 1.5 Mq (see e.g. m)- This implies that the GW emis¬ 
sion of the postmerger phase probes a density regime that 
is hardly accessible by other observational methods. In 
particular, attempts to infer NS properties via finite-size 
effects in the GW signal during the inspiral phase pre- 
ceeding the merger are restricted to the mass range of the 
individual components of the binary (see e.g. 157115^11^01 
I54U105U421I1261 I51 I501I55] '). In this sense the postmerger os- 

^ Within this paper the term “high density” refers to den¬ 
sities roughly above nuclear saturation density, which is the 
regime most relevant for the stellar structure of NSs. 


dilations represent a complementary approach by provid¬ 
ing insights to NS properties at higher masses and higher 
densities (e.g. [?101]). Moreover, the collapse behavior of 
NS merger remnants may reveal the properties of the very 
high density regime, which is decisive for determining the 
threshold of the gravitational collapse m- The consid¬ 
eration of the collapse behavior may be particularly re¬ 
warding for establishing limits on the maximum mass of 
nonrotating NSs, which is challenging given the paucity 
of NS systems with masses close to the collapse thresh¬ 
old, while the mass range of merger remnants partially 
overlaps with the mass range where the gravitational col¬ 
lapse is expected to take place for various EoSs. Therefore, 
an upper limit on the maximum mass may be established 
through NS mergers, which may be hard to derive through 
other types of observations. 

The basic impact of the EoS on the merger dynamics 
and GW signal can be described as follows. Stiff EoSs lead 
to NSs with large radii. Such stars can be deformed more 
easily under the influence of an external tidal field. Con¬ 
sequently, finite size effects during the inspiral set in at 
a larger orbital separation, i.e. a lower orbital frequency, 
and the stars finally merge at a relatively low orbital fre¬ 
quency. In contrast, soft EoSs yield more compact NSs, 
which in comparison to larger NSs behave more like point 
particles during the late inspiral phase. Stars described 
by soft EoSs are harder to deform and during the inspi¬ 
ral they reach higher orbital velocities before they merge. 
This also implies a higher linear velocity before merging 
and results in a merger with a higher impact velocity. The 
stiffness also affects the dynamics of the postmerger phase 
and, in particular, the frequencies of the excited oscilla¬ 
tion modes. A stiff EoS leads to a relatively large merger 
remnant, whose quadrupolar fluid oscillation frequency is 
relatively low, since it scales approximately with the mean 
density. In the case of a soft EoS, the merger remnant is 
more compact and thus oscillates at higher frequencies. 
(For such EoSs the higher impact velocity during merging 
additionally leads to a stronger excitation of the quasi- 
radial oscillation mode of the remnant.) 

In this paper we summarize the current status of possi¬ 
bilities to deduce NS properties and EoS constraints from 
GW observations of the postmerger phase of NS mergers. 
We focus mostly on the most prominent oscillation mode 
in the GW spectrum of the postmerger phase, the funda¬ 
mental quadrupolar fluid mode, since this is the feature 
which is most likely to be detected. We prepend a discus¬ 
sion of mass measurements through GW chirp signals of 
NS mergers, which is essential to NS radius measurements 
by the observation of GWs from the postmerger phase as 
well as for constraints on the maximum mass of nonrotat¬ 
ing NSs. We supplement a review of previous findings by 
new results, which for instance clarify the role of intrinsic 
NS rotation of the binary components. Also, we provide 
further evidence that the merger remnant can be consid¬ 
ered to be a single, self-gravitating object and can thus 
be investigated by concepts of asteroseismology By 
considering certain oscillation modes and secondary fea¬ 
tures of the postmerger GW spectrum, we arrive at a uni- 
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fied classification scheme of the postmerger dynamics and 
GW emission. We also present an analytic model of the 
GW emission, which may represent a step towards the 
construction of GW templates that can be employed in 
GW data analysis. Along the lines of the current Topical 
Issue of “Exotic Matter in Neutron Stars” we also dis¬ 
cuss a scenario of the existence of two families of compact 
stars [46| and identify certain observational features which 
can test these ideas. This discussion represents an exam¬ 
ple of how the postmerger phase can be employed to test 
high-density matter properties at densities which are not 
accessible through NSs of masses below ^ 1.5 Mq. 

The results presented in this paper rely mostly on 
the calculations discussed in [niiisiiiaiiHiiia, where fur¬ 
ther information can be found. Additional new findings 
laid out here are obtained within the same physical and 
numerical model, for which further details are provided 
in [MUMUmUKH]. We only consider mergers from quasi¬ 
circular orbits (for a discussion of tidal capture events see 
e.g. [82]). This paper is written in a modular way such that 
the individual sections are mostly self-contained and the 
reader can skip certain topics. Introductory remarks on 
mass measurements via GWs are presented in Sect. HI NS 
radius determinations via the postmerger phase are dis¬ 
cussed in Sect. |31 We address the collapse behavior of NS 
merger remnants and outline possibilities to measure the 
maximum mass of nonrotating NSs in Sect. 2] The origin 
of the major and of secondary features of the GW spec¬ 
trum are explained in Sect. [Sjas well as the dependencies 
of the secondary GW peaks, which lead to a classifica¬ 
tion scheme of the postmerger dynamics and GW emis¬ 
sion. In Sect. m we provide details on our analytic model 
of the postmerger GW emission. The impact of intrinsic 
NS rotation is considered in Sect.jTj The possible observa¬ 
tional implications of a scenario of two distinct families of 
compact stars is explored in Sect. |S| We summarize and 
conclude in Sect. M 


2 Mass measurements 

The strongest part of the GW signal of a NS merger origi¬ 
nates from the preceding inspiral phase, where the system 
continuously losses angular momentum and energy by the 
emission of GWs. The losses lead to inspiralling trajec¬ 
tories of the binary components, while the GW signal is 
determined by the orbital motion of the NSs. The inspiral 
accelerates, which results in a chirp-like GW signal with an 
increasing amplitude and an increasing frequency. Finite- 
size effects become important only for the last orbits be¬ 
fore merging and may yield information on the NS EoS, 
specifically on the tidal deformability, which for fixed mass 
approximately scales with the stellar radius (see e.g. [571 



The key parameter determining the GW signal during 
the inspiral phase is the so-called chirp mass apart from 
other parameters like the distance to the source, the in¬ 
clination of the binary, the orientation of the instruments 



Fig. 1 . Total binary mass Mtot = M1+M2 as a function of the 
mass ratio q = Mi/M 2 for systems with constant chirp mass 
(solid lines), which will be measured with high precision by GW 
detectors. Dashed lines show systems where the less massive 
component has a mass of Mi = 1.1 Mq and Mi = 1.2 Mq, 
which represent reasonable lower limits on the NS mass, and 
thus restrict the allowed systems to the right of the respective 
dashed curves. 

and so on. The chirp mass is given by 

Mehirp = + M2)-1/® (1) 

with the NS masses Mi and M 2 . Being the crucial quan¬ 
tity to determine the signal, the chirp mass will be mea¬ 
sured with very high accuracy since it is one of the param¬ 
eters used to construct the template bank for matched 
filtering searches that will be employed to actually de¬ 
tect GWs from binaries [551IMll751[TnilT^ITMll551[niTn71 
11251155] . In contrast, the mass ratio can only be measured 
with higher accuracy for nearby events |551l38l[73lll01ll231 
[T751l551[nil071ll251l55| . 

It can be seen from Eq. o that a measurement of 
the chirp mass alone does not suffice to determine the in¬ 
dividual binary masses, which can only be achieved by 
determining the mass ratio with sufficient precision. How¬ 
ever, for NS binaries the chirp mass alone provides already 
a very good estimate of the total binary mass Mtot = 
Ml -I-M 2 . Figure [1] shows the total binary mass as a func¬ 
tion of the mass ratio q = Mi/M 2 < 1 for sequences of 
constant chirp mass (solid lines). For a constant chirp mass 
the total binary mass depends only weakly on the mass 
ratio. The dashed lines define the binary systems, where 
the less massive component has a gravitational mass of 
1.1 Mq and 1.2 Mq, respectively. Under the reasonable 
assumption that NSs cannot be less massive than 1.2 Mq 
(see discussion below), a measured chirp mass restricts the 
possible binary parameters to the right of the dashed curve 
for Ml = 1.2 Mq. This implies that for instance a mea¬ 
sured chirp mass of 1.1752 Mtot restricts the total binary 
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mass of this event to the range 2.7 Mq < Mtot <2.73 Mq 
(without any information about the mass ratio). From ob¬ 
served NS binaries (see e.g. listing in [HI]) and theoretical 
population synthesis studies (e.g. gS]) it is expected that 
most binaries have a chirp mass close to 1.1752 Mtot (cor¬ 
responding to a total binary mass of roughly 2.7 Mq). 
This means that for the most likely merger event the to¬ 
tal binary mass can be determined very accurately from 
the chirp mass measurement only. Even for very massive 
binaries (with a chirp mass of 1.3058 Mq) the total binary 
mass can be inferred with a precision of ±0.05 Mq. 

To support the arguments above we remark that recent 
core-collapse simulations for a large number of progenitor 
stars find a minium NS rest mass of about 1.3 Mq ED- 
A NS rest mass of 1.3 Mq corresponds to a gravitational 
mass of roughly 1.2 Mq, somewhat dependent on the EoS. 
These findings justify to assume that NS binaries are un¬ 
likely to host NSs less massive than 1.2 Mq. Also other 
formation channels like an accretion-induced collapse of 
a white dwarf are unlikely to form less massive NSs. See 
also the discussion of the minimum NS mass in [81]. Ei- 
nally, we note that a white dwarf-NS binary might have a 
chirp mass in the range which is typical of NS-NS binaries, 
and a more extreme mass ratio if the white dwarf is less 
massive than 1.2 Mq. The merger of a white dwarf-NS 
binary will, however, lead to significantly altered inspiral 
dynamics m and thus a different GW signal prior to 
merging such that a white dwarf-NS merger can be easily 
distinguished from a NS binary merger. 

The considerations above play a role for the follow¬ 
ing discussions, because methods to infer EoS properties 
from NS mergers rely on the ability to measure at least 
the total binary mass from the inspiral GW signal. In¬ 
formation on the mass ratio and thus the individual bi¬ 
nary component masses is not critical, but improves the 
constraints on NS properties. As outlined above, the de¬ 
termination of the total binary mass via the chirp mass 
alone represents the absolute minimum of what is achiev¬ 
able with existing and upcoming GW detectors. However, 
it is important to stress that for distances of order 50 Mpc, 
which allow to deduce EoS properties, the individual bi¬ 
nary masses are expected to be recovered with a preci- 

I53| . For instance, in m the individual NS masses have 
been recovered within about 10 per cent at a distance of 
100 Mpc. Assuming that the error scales linearly with the 
distance, an accuracy of a few per cent could be reached in 
determining the individual masses of the merging NSs at 
a distance of a ^ 50 Mpc. Therefore, we will work in the 
following under the condition that the individual masses 
can be measured sufficiently accurately for our purposes 
and that, in particular, the total binary mass can be de¬ 
termined very well. 

3 Radius measurements 

The outcome of a NS merger essentially depends on the 
EoS of high-density matter and the total binary mass. The 
binary mass ratio has a modest impact, while the initial 



Fig. 2. GW spectrum of the cross polarization of a 1.35- 
1.35 Mq merger with the DD2 EoS |63II122J along the polar 
direction at a distance of 20 Mpc. heft = h{f) ■ f with the 
Fourier transform of the waveform hx and frequency /. /peak, 
/spiral and / 2-0 are particular features of the postmerger phase, 
which can be associated with certain dynamical effects in the 
remnant. Since the simulation started only a few orbits before 
merging, i.e. at a relatively high orbital frequency, the power 
at lower frequencies (below kHz) is massively underrepre¬ 
sented in the shown spectrum, and the low-frequency part of 
the spectrum does not show the theoretically expected power- 
law decay. The thin solid lines display the spectra of the GW 
signal of the postmerger phase only revealing that the peaks 
are indeed generated in the postmerger phase. Dashed lines 
show the expected unity SNR sensitivity curves of Advanced 
LIGO |61| (red) and of the Einstein Telescope |66 | (black). 


NS rotation and magnetic fields have only a small influ¬ 
ence.) Simulations show that for binary masses of about 
2.7 Mq the merging leads to the formation of a massive, 
hot, differentially rotating NS remnant for most EoSs, 
even for EoS models which yield a maximum mass of non¬ 
rotating NSs significantly below 2.7 Mq, e.g. |130lll09l[m 
i?mii 11 1ll 1 3111 1 2i[ieiii?im[Mirr^[7miTT^i^[7niiTTTniT7ir^ 

[T7MIT^[Ml[^l^[TT7l[7^ITTMITT^l74l[^IMllTTmi^ 

The rapid rotation stabilizes the remnant against gravi¬ 
tational collapse. According to pulsar observations in NS 
binaries most systems are expected to have a total mass of 
roughly 2.7 Mq (see e.g. the compilation of binary masses 
in [81]), which is also supported by theoretical studies of 
the binary population, e.g. |45] . Given this mass range, 
the most likely outcome of a merger event should be ex¬ 
pected to be the formation of a NS remnant (see [15] for 
the EoS dependence of collapse behavior). 

A typical GW spectrum of a NS merger resulting in 
the formation of a NS remnant is shown in Eig. [5] The 
spectrum is computed from a simulation of a 1.35-1.35 Mq 
merger with the DD2 EoS [5511122] . The dominant GW os¬ 
cillation frequency /peak of the postmerger phase is clearly 
visible as a pronounced peak in the kHz range, which is 
produced by the dominant quadrupolar remnant oscilla¬ 
tion. Apart from the main peak one recognizes several 
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additional peaks, whose nature will be discussed in Sec¬ 
tion By computing the spectrum of the GW signal of 
the postmerger phase only, one can show that the differ¬ 
ent features are related to the postmerger phase (see thin 
lines in Fig. [5]). 

The dominant oscillation frequency /peak depends in 
a specific way on the high-density EoS [TfllTB] . This is 
understandable since the EoS affects the size of the rem¬ 
nant, which in turn determines its oscillation frequency 
(see Fig. 13 in [15]). The structure of the remnant is 
also influenced by its angular momentum. The available 
angular momentum, however, is given by the dynamics 
of the late inspiral/merging phase, which is fully deter¬ 
mined by the stellar structure of the inspiralling stars and 
thus also depends on the EoS in a particular way. The 
strong EoS dependence of the peak frequency can be ex¬ 
pressed as follows. An EoS which is used in a given sim¬ 
ulation, can be conveniently characterized by the radii 
of nonrotating NSs, which are uniquely determined by 
this EoS through the stellar structure equations (Tolman- 
Oppenheimer-Volkoff equations |121l[M] l. Specifically, for 
a set of calculations with a fixed total binary mass but 
different EoSfl, we relate the peak frequency, which is 
extracted from a simulation with a given EoS, to the ra¬ 
dius of a nonrotating NS (described by the same EoS) 
with a fixed fiducial mass. A natural choice is to employ 
the NS radius for a mass of Mns = Aftot/2, which for 
symmetric binaries is just the radius of the inspiralling 
NSs (more precisely, at infinite orbital separation). In this 
case a clear correlation is found, where EoSs leading to 
more compact NSs yield higher postmerger GW frequen¬ 
cies (see Eig. 12 in [TB], which shows this relation for 
Mtot = 2.7 Mq). (Alternatively, one can use the com¬ 
pactness C = GMNs/(c^.R(AfNs)) of fiducial nonrotating 
NS models, which is equivalent to employing the radius 
i?(MNs)-) 

The empirical relation between /peak and R{M^^) is 
very tight, which implies that a measurement of the peak 
frequency can be used to determine the unknown radius 
of a nonrotating NS with a fixed mass by simply inverting 
the empirical relation HZEB]. Thus, a future detection of 
the GW postmerger phase and extraction of the peak fre¬ 
quency (see [551136] ! will yield strong constraints on the 
high-density EoSs. In [TTIITS] the largest deviation of the 
empirical data from a fit is only a few hundred meters. 
The accuracy of a radius determination by the postmerger 
GW signal is mostly affected by two sources of error. One 
error is the uncertainty of the measurement of the peak 
frequency. Apart from this, one should take into account 
deviations between the data and the fit to the data de¬ 
scribing the empirical relation. A measurement of the peak 
frequency (of the true EoS) does not reveal in which way 


^ Except for some models considered in Sect. 21 the EoSs 
discussed in this study are temperature dependent and include 
electrons, positrons and photons, while neutrino contributions 
are neglected. With regard to the resulting stellar properties 
these EoSs cover a representative range, which, for example, 
can be seen from the range of radii in Fig. [3] and maximum 
masses in Fig. [8] 



Fig. 3. Dominant postmerger GW frequency /peak as function 
of the radius Ai.e of a nonrotating NS with a gravitational mass 
of 1.6 Mq for different EoSs and different total binary masses 
(plus signs for 2.4 Mq, circles for 2.7 Mq, crosses for 3.0 Mq) 
and a mass ratio of unity. The solid lines are least-square fits 
to the data of the different binary masses. 



Fig. 4. Dominant postmerger GW frequency /peak as function 
of the radius Ai.e of a nonrotating NS with a gravitational mass 
of 1.6 Mq for symmetric and asymmetric binaries with a chirp 
mass of Mchirp = 1.1665 Mq for different EoSs. Black symbols 
display data for 1.34-1.34 Mq mergers, while blue symbols ex¬ 
hibit peak frequencies of 1.2-1.5 Mq mergers. The solid line is 
a least-square fit to the data of the different binary masses. 


the measured frequency slightly deviates from the empiri¬ 
cal relation. Hence, one conservatively has to assume that 
the true data point may deviate as much as the largest 
deviation found in the large sample of candidate EoSs. 

The peak frequency has been shown to be measurable 
with very high precision by a coherent burst search anal¬ 
ysis |35j . In this study waveforms from numerical models 
were superimposed with the recorded data stream of pre¬ 
vious GW detector science runs, which simulates the noise 
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of the future instruments. The model waveforms were in¬ 
jected at random times and the noise was rescaled to the 
anticipated sensitivity of the second-generation GW detec¬ 
tors Advanced LIGO and Virgo. The existing GW data 
analysis pipeline was able to recover the injected signal 
and to determine the peak frequency with an accuracy of 
^10 Hz, which is smaller than the spread in the empiri¬ 
cal relation between /peak and the NS radius. This implies 
that the radii of the inspiralling stars can be determined 
with a precision of a few hundred meters. 

These considerations show that the larger contribu¬ 
tion to the error of a radius measurement originates from 
the scatter in the empirical relation between /peak and 
R{Myis)- In this context, the following observation is im¬ 
portant. One has the freedom to choose any fiducial NS 
mass different from Mns = Vftot/2 for characterizing a 
given EoS by the TOV radius i?(MNs)- Empirically, it 
turns out that using a fiducial NS mass somewhat larger 
than Mns = Mtot/2 leads to tighter relations between 
/peak and i?(MNs)- This is exemplified in Fig.|31 For Mtot = 
2.7 Mq (circles in Fig. [3]) the maximum deviation be¬ 
tween the data and a fit amounts to only ^ 175 me¬ 
ters if Mns = 1-6 Mq is chosen. This implies that the 
measurement of the dominant postmerger frequency for 
Mtot = 2.7 Mq determines the radius of a nonrotating 
1.6 Mq NS with an accuracy of better than 200 meters. 

It is natural that a fiducial mass of Mns = 1-6 Mq 
is somewhat more appropriate than Mns = 1-35 Mq for 
characterizing the postmerger oscillations of 1.35-1.35 Mq 
mergers ( Mtot = 2.7 Mq). The maximum densities in the 
massive, rotating merger remnant are higher than in the 
initial NSs and they are comparable to the central den¬ 
sities of nonrotating, static NSs with a mass of roughly 
1.6 Mq (see e.g. Fig. 15 in [l^)- For this reason, nonro¬ 
tating NSs with Mns > Mtot/2 better represent the den¬ 
sity regime encountered in the merger remnant and thus 
provide a better description of the EoS. 

Figure [3] also shows that similar empirical relations 
hold for other binary masses. The different symbols dis¬ 
play the peak frequencies from calculations with total bi¬ 
nary masses of 2.4 Mq (plus symbol), 2.7 Mq (circles), 
and 3.0 Mq (crosses) for symmetric mergers. The solid 
lines are least-square fits to the data. One recognizes that 
also for other binary masses there are only small devi¬ 
ations from the empirical relations on the order of only 
a few 100 meters. Hence, these relations can be used for 
determining NS radii with similar precision after the to¬ 
tal binary mass was deduced from the inspiral signal (see 
discussion in Sect. HI. The total binary mass determines 
which of the relations has to be employed to convert the 
measured peak frequency to a NS radius. Note that all 
data are plotted as a function of i?(1.6 Mq) and hence 
the fiducial mass is not optimized to yield the tightest re¬ 
lations for total binary masses of 2.4 Mq or 3.0 Mq. Here 
we assume that the determination of the binary masses 
does not contribute a significant error, but that the bi¬ 
nary masses are known with sufficient precision as argued 
in Sect. [H We also point out that very similar relations 


can be constructed for unequal-mass binaries with a fixed 
total mass and mass ratio. 

We here do not further discuss asymmetric binaries in 
great detail, but refer for instance to [18]. There it is shown 
that the peak frequencies of asymmetric binaries deviate 
only somewhat from the ones of the symmetric binaries 
of the same total mass. This can be also seen in Fig. [4] 
which shows the peak frequency as function of Ri,e for 
symmetric and asymmetric mergers. The symmetric and 
asymmetric binary systems are chosen such that they have 
the same chirp mass of Mchirp = 1.1665 Mq (see Sect.[3|), 
which will be measured very accurately from the GW in¬ 
spiral signal. This chirp mass corresponds to a symmetric 
setup with two stars of 1.34 Mq (black symbol^ in Fig. 11, 
i.e. Mtot = 2.68 Mq (see Eq. CD). The asymmetric config¬ 
urations with Mchirp = 1.1665 Mq shown as blue symbols 
in Fig. m are 1.2-1.5 Mq mergers, which are arguable the 
most asymmetric systems expected in the binary popula¬ 
tion (see discussion in Sect. HD- Less asymmetric binaries 
lead to smaller deviations from the peak frequency of the 
symmetric system. Importantly, the data in Fig. |4| illus¬ 
trate that even for the unlikely case that no information 
on the binary mass ratio is available and only the chirp 
mass was measured (Sect. [3|), the dominant postmerger 
frequency /peak still determines NS radii very accurately. 
The data points deviate by at most 258 meters from a 
fit (solid line), which is constructed from the data of the 
symmetric and asymmetric binaries. 

The precise form of the fits of the /peak — R data in 
Fig. [3| and the deviations from the fits of course depend 
slightly on the chosen sample of candidate EoSs, which 
should preferentially include all EoSs that are acceptable 
based on current knowledge. Moreover, uncertainties in 
the numerical or physical model may also lead to slight 
differences in the empirical relations, which, however, will 
diminish in the future by more detailed simulations. Also, 
it is possible to choose different functional forms which 
may to some small extent affect the accuracy of the empir¬ 
ical relations. It is worth noting that in contrast to m, 
we do not see difficulties in fitting a function to their data 
for the dominant oscillation frequency. Using the listed 
frequencies and computing the radii of 1.6 Mq NSs for 
the employed piecewise polytropic EoSs, it is easily possi¬ 
ble to fit for instance a bi-linear (or quadratic) function to 
the data for a fixed Mtot- This, for instance, results in de¬ 
viations of at most 214 meters for the 1.3-1.3 Mq mergers, 
and thus confirms our conclusions that a measurement of 
/peak is sufficient for a radius determination. 

In this context we note that the peak frequencies from 
our numerical simulations agree well with those of |lini 
13311 11 711 11 81153] . Hence, it would be surprising if the data 
of |117U118j did not yield a close relation for a fixed bi- 

® We did not simulate 1.34-1.34 Mq mergers, but interpo¬ 
lated linearly the /peak(Mtot) relation, which is given by the 
results from 1.2-1.2 Mq binaries and 1.35-1.35 Mq binaries. 
In this binary mass range a linear interpolation is a very good 
approximation for symmetric systems (see Fig. 1 in |23p. 

■* Note that the dominant oscillation frequency is called /2 
in as well as in some other references. 
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Fig. 5. Rescaled dominant postmerger GW frequency 
/peak/A/tot as function of the radius Ri.e of a nonrotating NS 
with a gravitational mass of 1.6 Mq for different EoSs and 
different total binary mass (plus signs for 2.4 Mq, circles for 
2.7 Mq, crosses for 3.0 Mq) and a mass ratio of unity. 

nary mass while our data show a good correlation and the 
frequencies essentially agree. 

In [TB] it was argued that relations between /peak and 
NS radii are not unexpected because the dominant emis¬ 
sion is generated by the fundamental quadrupolar oscilla¬ 
tion mode m , whose frequency is known to scale tightly 
with ^JM/R^ for nonrotating NSs (see 0)- This scaling 
for nonrotating stars suggests that it may be possible to 
rescale the frequency by dividing by i/Mtot to describe 
the data of different total binary masses by a single re¬ 
lation. While this is indeed possible, empirically we find 
that fpeak/Mtot results in an even tighter relation. Fig¬ 
ure [S] shows the rescaled peak frequency /peak/Aftot as a 
function of Ri,e = 77(1.6 Mq) for total binary masses of 
2.4 Mq, 2.7 Mq, and 3.0 Mq. The rescaled frequencies 
can be described by the quadradic least-squares fit 

/peak/Mtot = 0.0157 • Rle - 0.5495 • Ri.e + 5.5030, (2) 

with /peak in kHz, the binary mass in Mq and 7?i.6 in km. 
The data in Fig. deviate by at most 500 meters from 
the fit. The reason why fpeak/Mtot yields a better uni¬ 
versality than /peak /VMtot (which results in deviations of 
up to ^ 1 km) may be that the radius is also mass de¬ 
pendent. A perfect scaling behavior with the mass cannot 
be expected considering the findings in [23], which show 
that different EoSs can lead to different dependencies on 
the total binary mass (see Fig. 1 in [13]) • While such a 
universal relation of the rescaled /peak is theoretically in¬ 
teresting, from a practical point of view it is preferable to 
consider different sets of mass-dependent relations as in 
Fig. 131 (or possible interpolations between them) because 
they yield tighter relations, and the total binary mass will 
always be known from the GW inspiral signal. Concerning 
mass-independent relations for /peak we also note that a 
relation between the dominant postmerger frequency and 


the tidal coupling constant, pointed out in m , is inter¬ 
esting because it connects the inspiral and the postmerger 
signals. Such a relation is understandable because essen¬ 
tially both quantities are known to depend on the EoS 
and to scale with the NS radius. As in the case introduced 
above (Fig. |S]) the practical use of such relations depends 
on their tightness and a spread of ^ 500 Hz may be large 
compared to a typical width of the main postmerger peak. 

For a first assessment of the detectability of the domi¬ 
nant postmerger GW frequency by a morphology-independent 
burst search data analysis, see [33]. A more detailed study 
employing a principal component analysis and evaluating 
the detectability for discussed future detectors is presented 
in [33]. 


4 Gravitational collapse and estimates of the 
maximum mass 

Apart from determining NS radii as discussed in the previ¬ 
ous section, the postmerger phase also offers the possibility 
to constrain the maximum mass of nonrotating NSs. The 
maximum mass is a key parameter regarding the proper¬ 
ties of high-density matter because its value is determined 
by the very high-density regime that is usually not encoun¬ 
tered in most observed NSs, which have a lower mass. The 
maximum mass is thus particularly interesting for prob¬ 
ing exotic phases of matter, which potentially occur at 
higher densities. Here we outline three different methods 
to determine or at least constrain the maximum mass of 
nonrotating NSs. An important quantity in the following 
discussion is the threshold binary mass to prompt black 
hole (BH) collapse. If one considers binary mergers with 
different total binary masses for a given EoS, one finds 
that for a binary mass above a certain threshold the merg¬ 
ing leads directly to the formation of a BH and no (tran¬ 
siently stable) NS remnant is formed because the merged 
object is unstable against gravitational collapse. Note that 
the gravitational collapse in a merger remnant sets in at 
densities below the maximum density of non-rotating NSs. 

The connection between this threshold mass and the 
maximum mass of nonrotating NSs has the following back¬ 
ground. The maximum NS mass M^ax is the threshold 
for BH formation of static, nonrotating stars, while the 
threshold mass Mthres represents the threshold for BH for¬ 
mation of differentially rotating, hot NSs. Thus, a quan¬ 
titive relation between those two quantities is expected; 
and a determination of Mmax is possible since the thresh¬ 
old mass can be determined from the binary masses of 
different merger events considering their respective out¬ 
come. 

Method 1: An analysis of a large number of simulations 
with different EoSs and different total binary masses has 
revealed that the threshold binary mass Mthres depends in 
a particular way on the EoS [15] . This dependence can be 
described by TOV properties, which are uniquely defined 
by the EoS. With very good accuracy the threshold mass 
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Fig. 6. Threshold binary mass for prompt BH formation as 
function of the maximum mass of nonrotating NSs for different 
Ri.6- The solid lines are given by Eq. which describes the 
empirical data with a good accuracy m- 



2.6 2.8 3 


3.2 3.4 


Fig. 8. Maximum mass of nonrotating NSs as function of 
the dominant postmerger GW frequency /peak for 1.5-1.5 Mq 
mergers with different EoSs. The blue symbols mark EoSs 
which are particularly stiff at higher densities. The dashed line 
roughly indicates an upper limit on Mmax for a given peak 
frequency, see Eq. ©• 



Fig. 7. Threshold binary mass for prompt BH formation as 
function of the dominant postmerger GW frequency /peak of 
1.5-1.5 Mq mergers for different EoSs. 


is given by 

Mthres = (^-3.606 ^^^^"^" + 2.38^ • (3) 

with the gravitational constant G and the speed of light 
c. Similarly, it can also be expressed as 

Mthres = (-3.38^^ + 2 . 43 ) • M^ax, (4) 

i.e. using the maximum mass and the compactness of the 
maximum-mass configuration. The coefficients in both equa¬ 
tions are obtained by fits to the ratio , where Mthres 


is given by the results from simulations m- Equations (| 31 ) 
and reproduce the numerical results with a precision 
better than 0.1 Mq (on average the deviations between the 
numerical and the estimated threshold mass are 0.024 Mq 
for Eq. ([3]) and 0.031 Mq for Eq. (|4])). Bear in mind that 
the numerical value of Mthres can only be determined with 
a finite precision because apart from uncertainties of the 
numerical and physical model only a limited set of mod¬ 
els with different Mtot has been analyzed. m tested the 
collapse behavior for two EoS models, which were also 
employed in [15], and found compatible results. 

Using Eq. o Fig. |6| visualizes how Mthres depends on 
the maximum mass of nonrotating NSs for different cho¬ 
sen values of If the radius Ri,e has been determined 
observationally, for instance through the detection of the 
peak frequency of a 1.35-1.35 Mq merger (see Sect. [3|), 
then the threshold mass becomes an unambiguous func¬ 
tion of Mmax- The threshold mass may be observation- 
ally determined by confirming or excluding the presence 
of postmerger GW emission of a NS remnant for several 
merger events with different total binary masses. (Note 
that the prompt formation of a BH leads to very weak 
postmerger GW emission, which can be distinguished from 
the formation of a NS remnant for near-by events and suf¬ 
ficient sensitivity.) The inversion of the relation between 
Mthres and M^ax for a fixed Ris can then be employed 
for an estimate of the maximum mass of nonrotating NSs. 

We remark that a distinction between the prompt col¬ 
lapse and the formation of a NS remnant may also be 
possible by the observation of electromagnetic counter¬ 
parts of mergers. In particular, the thermal emission of 
the ejecta, which are heated by radioactive decays, may 
be observable The ejecta mass depends on 

the EoS, binary mass and mass ratio (see e.g. [in])- The 
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prompt collapse to a BH leads to smaller amounts of un¬ 
bound matter compared to mergers which form a NS rem¬ 
nant. This should lead to a noticeable difference in the 
light curves of these two different outcomes and may allow 
an observational identification of prompt collapse events. 
The detectability of electromagnetic counterparts to GW 
events has been addressed for instance in [QOllMlIS^ and 
depends on (currently) uncertain details like the observer 
angle, the available instruments, the search strategy, the 
observing conditions, the ejecta mass, the outflow veloc¬ 
ity, and the neutron-richness of the ejecta determining the 
opacity. The last three source properties are affected by 
the high-density EoS and the binary parameters. 

Method 2: We point out that already a single GW event 
may provide a relatively precise estimate of the threshold 
mass and a constraint on the maximum mass of nonrotat¬ 
ing NSs. Figure [7] displays the threshold mass as a func¬ 
tion of the peak frequency of equal-mass mergers with 
Mtot = 3.0 Mq for different EoSs. We use here the thresh¬ 
old masses estimated via Eq. which are very close 
to the numerical values. The blue symbols show the re¬ 
sults for additional simulations with representative EoSs 
from [62) . which are derived within a chiral effective field 
theory at lower densities (up to roughly nuclear saturation 
density) and supplemented by extrapolations of the EoS 
at higher densities with physically motivated parameter 
variations (minimum NS mass of 2 Mq and causality). 
For more information on this set of EoSs see [TH]. Since 
these EoSs describe only cold NS matter, the EoSs are 
supplemented with an approximate treatment of thermal 
effects for the merger simulations (see e.g. [H] for details 
and an assessment of this approximation). In these calcu¬ 
lations the “ideal-gas” index Tth, which regulates the ther¬ 
mal pressure support, was chosen to be 5/3. See e.g. [TO] 
and m for a motivation of this value. One should keep 
in mind that the exact choice of Tth introduces an am¬ 
biguity and thus an uncertainty of the numerical values 
of certain quantities extracted from simulations. This is 
the reason why we did not include these models in the 
discussion in Sect. [31 Note that from the original set of 
six EoSs in m. two models lead to the prompt formation 
of a BH for 1.5-1.5 Mq binaries. This is fully compatible 
with the theoretically estimated threshold mass from the 
TOV properties of these EoSs via Eqs. © and 0 . 

Figure [7] shows that the determination of /peak of a 
1.5-1.5 Mq merger yields the threshold mass with a pre¬ 
cision of about 0.2 Mq. Note that a total binary mass of 
3.0 Mq is only slightly larger than the masses of observed 
binaries for which precise mass measurements exist. Thus, 
a detection of a merger with Mtot = 3.0 Mq may not be 
improbable. 

At first glance it may seem tempting to use the peak 
frequency of a 1.5-1.5 Mq merger to fix Mthres (Fig. [7]), 
then to determine Ri,e as discussed in Sect. [3] (Fig. [3]), and 
then to estimate iWmax via Eq. This, however, repre¬ 
sents an attempt to determine two independent quantities 
{Mmax and Ri.e) by only one observable. In fact, a thor¬ 
ough error analysis for only one peak frequency measure¬ 
ment reveals that depending on the EoS the maximum 


mass can be determined only with a large error bar. (No¬ 
tice the flat slope for small radii in Fig. 1) These con¬ 
siderations are summarized in Fig. [51 which shows Mmax 
directly as a function of the peak frequency of 1.5-1.5 Mq 
mergers. Remarkably, /peak(Altot = 3.0 Mq) immediately 
yields an upper limit for Mmax, which for soft EoSs implies 
even a good estimate of M^ax (recall the lower bound on 
Afmax by the measurement of NSs with about 2 Mq). In 
Fig. m the dashed line, which is given by 


A/jnax.upperlimit — ^fpea]i T 4.53, 


(5) 


with the frequency in kHz and the mass in Mq, indicates 
an exclusion region and thus defines an upper limit on 
Afmax, which can be derived from a single detection of a 
merger with Mtot = 3.0 Mq. We stress that some of the 
EoSs within our sample (especially some of the models 
by [SI], blue symbols) have a speed of sound equal to 
the speed of light at higher densities. These models are 
thus maximally stiff, which is favorable for yielding high 
maximum masses. We thus conjecture that the dashed line 
represents a true upper limit and is not an artifact of the 
chosen sample of EoSs. We also refer to Fig. 16 in |18) . 
which shows a plot similar to Fig. jSjbut for simulations of 
1.35-1.35 Mq mergers. 

Fig. |5| shows that a single merger event can only con¬ 
strain the maximum mass. However, we emphasize that 
the estimate of the threshold mass, e.g. through Fig. (TJ 
represents a viable goal on its own. The threshold mass 
is the crucial quantity to judge the outcome of a detected 
merger event. For instance, at larger distances only the 
chirp mass can be measured (which determines the total 
mass with some accuracy, see Sect. la. Having an estimate 
of Mthres available from a previous near-by event, the to¬ 
tal mass of a distant merger may be crucial information 
for interpreting a possibly detected electromagnetic coun¬ 
terpart, for instance from thermal emission by the ejecta, 
which is powered by the radioactive decay of nucleosyn¬ 
thesis products |831l79l[9ni89] . It is known that the prompt 
collapse leads to a reduced ejecta mass and thus to dim¬ 
mer electromagnetic counterparts, which reach their peak 
emission on a shorter time scale e.g. m- Also, for the co¬ 
incident detection of a short gamma-ray burst [MiisniiiT] 
and a GW signal of a merger the information on the col¬ 
lapse behavior will provide valuable information and help 
to understand the conditions leading to this gamma-ray 
burst. 

Method 3: While a single event can only yield an up¬ 
per limit on the maximum mass M^ax, we have recently 
shown that several events with a measurement of the dom¬ 
inant postmerger GW frequency may suffice to determine 
the maximum mass m- Specifically, mergers with differ¬ 
ent binary masses (but in the most likely range of binary 
masses) can be employed for an M^ax estimate. The pro¬ 
cedure has the following background. Measuring the peak 
frequencies for two different total binary masses allows 
to estimate how the dominant GW frequency depends on 
Mtot- This can be used to extrapolate /peak(Mtot) and 
thus to estimate the behavior at higher binary masses, 
which is particularly sensitive to the maximum mass and 
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the radius of the maximum-mass configuration. In essence, 
the absolute values of /peak(-^tot) and its slope are very 
sensitive to M^ax and i?max and thus yield an estimate 
of the maximum mass of nonrotating NSs with an ac- 
cuarcy of roughly ±0.1 Mq. For details we refer to [53] and 
note that a slightly different approach is described in [53] , 
which results in a similar precision. We also point out that 
this extrapolation procedure determines the radius of the 
maximum-mass configuration and the maximum density 
of nonrotating NSs, which are quantities that are highly 
characteristic for the high-density regime of the EoS. 

The advantage of this approach lies in the fact that it 
relies only on detections of mergers with binary masses in 
the most likely range, i.e. roughly between 2.4 Mq and 
3.0 Mq (see e.g. [SI])- This contrasts the procedure to 
estimate iWmax by a direct measurement of the threshold 
mass as sketched for Method 1. Estimating the thresh¬ 
old mass directly by determining the outcome of mergers 
with binary masses above and below Mthres requires sev¬ 
eral merger detections at higher binary masses, which are 
possibly very unlikely. 

5 Characterization of GW peak frequencies 

5.1 Origin of main and secondary GW peaks 

We start with a brief summary of our current under¬ 
standing of the mechanisms which produce certain peaks 
in the GW spectrum For the interpretation of the 
main and secondary peaks it is important to realize that 
the merger remnant can be considered as an isolated self- 
gravitating object mg. This view is supported by Fig. [9] 
showing the lapse function in the equatorial plane for a 
1.35-1.35 Mq merger with the DD2 EoS [H31I122] . The 
lapse function may be seen as the relativistic analog of 
the gravitational potential. In the upper left panel the 
system is still composed of two cores corresponding to the 
initial NSs. Shortly later, the two cores have merged into 
a single potential well. The time steps of the snapshots in 
Fig. [g are the same as the ones in Fig. 2 of m, which 
shows the evolution of the rest-mass density in the equa¬ 
torial plane for the same model. It is remarkable that in 
the density evolution, a double-core structure persists for 
several milliseconds after merging, whereas the lapse func¬ 
tion exhibits already a single core with partially strong 
deformations. (In the density plots corresponding to the 
upper right and lower panels of Fig. j^a clear double-core 
structure is still visible.) These findings imply that the 
double-core structure visible in the density field for many 
milliseconds should not be interpreted as two gravitation¬ 
ally interacting objects, but rather as local overdensities 
moving in a single gravitational potential. Hence, the dou¬ 
ble cores should be viewed as tracers of the dynamics of 
a single, isolated object, rather than two independent dy¬ 
namical objects. 

The time evolution of the remnant structure is also 
summarized in Fig. 1101 It shows the time evolution of the 
central lapse function, which may be interpreted as a mea¬ 
sure for the compactness of the stellar object. The two ini¬ 


tial NSs first touch at about t = 12.4 ms. The vertical line 
indicates the time when the two cores in the lapse func¬ 
tion (upper left panel in Fig. [g merge into a single core, 
which occurs already during the first compression phase, 
i.e. right during the final plunge. The subsequent oscil¬ 
lations in the central lapse function indicate quasi-radial 
oscillations (bounces and compressions of the remnant). 
The figure implies that matter accumulates into a sin¬ 
gle, gravitational trough already at a very early time in 
the remnant’s evolution. Consequently, right from its for¬ 
mation the remnant’s evolution can be described by the 
(non-linear) dynamics of a single, self-gravitating object. 
For instance, the dominant osillation can be associated 
with the fundamental quadrupolar fluid oscillation mode. 
This can be seen by adding a velocity perturbation to the 
remnant at late times, when it reaches a quasi-equilibrium 
and the GW emission essentially has diminished. We add 
an instantaneous perturbation of 

Svg = 0.4sin (TT- ) sin (0) cos (0) cos (2(^), (6) 

\ ^surface (p) / 

to the 0-component of the coordinate velocity as a func¬ 
tion the polar angle 0, the azimuthal angle (j) and the radial 
coordinate r. (Geometrical units are adopted.) The radial 
coordinate Tsurface of the surface, which only depends on 
0, is defined by the coordinate at which the density drops 
below some threshold. We evolve the remnant with the 
added velocity perturbation for several milliseconds and 
extract the GW signal. The spectrum of the GW signal of 
the orginal simulation is compared to the spectrum of the 
perturbed model in Fig. [TT] It is evident that the pertur¬ 
bation excites an oscillation mode with the same frequency 
as the dominant remnant oscillation, which strongly sug¬ 
gests that /p^k is the frequency of the / = |to| = 2 funda¬ 
mental modq3. This is corroborated by the extraction of 
the oscillation eigenfunction in [ng, which shows a clean 
quadrupolar structure. 

Also other (secondary) peaks in the GW spectrum can 
be explained by oscillation modes of the remnant. For in¬ 
stance, the peak at 1.5 kHz in Fig. [5] is produced by a 
quasi-linear interaction between the fundamental quadrupo¬ 
lar mode and the quasi-radial oscillation of the remnant mg. 
This can be shown by determining the frequency /o of 
the quasi-radial mode, which itself does not occur in the 
GW spectrum. However, it is very pronounced in the time 
evolution of the central lapse function (Fig. [ig) and in 
other characteristic properties of the remnant, such as the 
maximum density and the size of the remnant. Again, by 
adding a suitable velocity perturbation to the late-time 
remnant, one can predominantly excite the quasi-radial 
mode and can extract its frequency from the time evolu¬ 
tion of the central lapse function. This is visible in Fig.lTUl 
where the green dashed curve shows the evolution of the 
radially perturbed model. The determination of the fre¬ 
quencies of the quadrupolar mode (/peak) and of the ra¬ 
dial mode (/o) reveals that a secondary peak is expected 

® For this reason it was called /2 in mg in the context of 
mode identification. 
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Fig. 9. Evolution of the lapse function of a 1.35-1.35 Mq merger with the DD2 EoS in the equatorial plane. The snapshots 
should be compared to Fig. 2 in showing the density evolution for the same time steps of the same simulation. 


to occur at /peak — /o) which indeed is the case. The fre¬ 
quency coincidence confirms the nature of this secondary 
peak as being a coupling of two modes, which is why we 
refer to this feature as the / 2-0 peak. Another correspond¬ 
ing combination frequency can be recognized in Fig. [5] at 
approximately /peak + /o- The peak at /peak + fo is obser- 
vationally less interesting because of its weakness and the 
smaller sensitivity of GW detectors at higher frequencies, 
but it substantiates the importance of mode couplings in 
NS merger remnants. 

Finally, there is one more secondary peak visible in the 
GW spectrum displayed in Fig. [5] Recently, we provided 
evidence that this feature is generated by a spiral defor¬ 
mation which is created during merging. This deformation 
cannot follow the faster rotation of the inner remnant. 
The spiral deformation forms antipodal bulges, which or¬ 
bit around the central part of the remnant for several mil¬ 
liseconds. Being a strong, non-axisymmetric, orbiting de¬ 
formation, the antipodal bulges generate a GW signal at 
a frequency which is twice the orbital frequency of the 


bulges. A deeper analysis of the simulation data confirms 
this origin of the GW peak at 2 kHz for teh particular 
model (e.g. by extracting the orbital motion of the an¬ 
tipodal bulges and estimating their mass, by comparing 
the presence of the bulges and the presence of the sec¬ 
ondary peak in the GW spectrum, and by computing GW 
spectra for the inner and outer parts of the remnant sep¬ 
arately to estimate the contribution of the different rem¬ 
nant components to the different GW features). See [35] 
for more details. Clearly, this so-called /spiral feature can¬ 
not be explained within a perturbative approach. One can 
show by an analytic model that a peak with an appropri¬ 
ate strength in the GW spectrum can be produced by 
point particles of a few 0.1 Mq, which orbit for only a 
few milliseconds with an orbital separation of roughly the 
diameter of the inner remnant (see also Sect. (5]). 
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Fig. 10. Time evolution of the central lapse function of a 1.35- 
1.35 Mq merger with the DD2 EoS (black line). The vertical 
line marks the time from which on only a single minimum of 
the lapse function is present. The dashed green line shows the 
evolution of the central lapse function of a late-time remnant 
of the same model where a quasi-radial velocity perturbation 
was added to excite the fundamental quasi-radial oscillation 
mode. 



Fig. 11. GW spectrum of a 1.35-1.35 Mq merger with the 
DD2 EoS (black line) given by hefi = h{f) ■ f with the Fourier 
transform of the waveform hx- The green dashed curve shows 
the GW spectrum of a simulation of a late-time merger rem¬ 
nant of the same model which was perturbed with a velocity 
field suitable to excite the fundamental quadrupolar fluid os¬ 
cillation mode. Thin dashed lines show the anticipated unity 
SNR sensitivity curves of Advanced LIGO [61] (red) and of the 
Einstein Telescope [SS] (black). 


5.2 Classification of postmerger GW emission 

Considering models of NS mergers with varied binary masses 
and with different EoSs, one realizes that the three most 
prominent features in the GW spectrum can be explained 
by the three mechanisms detailed above: the fundamental 
quadrupolar mode, the coupling of the quadrupolar and 
the quasi-radial mode, and the orbital motion of antipodal 
bulges. The peak of the fundamental mode is present in 
all models and it is always the strongest feature. The pres¬ 
ence and the strength of the different secondary features, 
however, are sensitively affected by the total binary mass. 
Depending on the total binary mass relative to the thresh¬ 
old mass Mthres (see Sect. n the different secondary peaks 
are more or less pronounced. This leads to a classification 
scheme that relies on the presence and strength of the dif¬ 
ferent secondary features. One can identify three different 
types of NS mergers and corresponding GW spectra. 

For relatively high Mtot, i-e. close to but below Mthres, 
the quasi-radial mode is strongly excited during merging, 
and consequently the / 2-0 feature is the most prominent 
secondary feature (Type I). In comparison, the /spiral peak 
is weaker and may even be hardly visible in the spectrum. 
For moderate binary masses both secondary features are 
clearly visible and can be clearly distinguished in the spec¬ 
trum (Type II). The secondary peaks have a roughly com¬ 
parable strength. In a third type of GW spectra, the / 2-0 
feature is absent (or hardly visible) and the GW peak by 
the spiral deformation is the most prominent secondary 
feature. This Type III occurs for relatively low total bi¬ 
nary masses, that means for Mtot much below Mthres- (An 
additional case is the prompt collapse to a BH as discussed 
in Sect. mfor Mtot > Mthres-) 

Since the threshold mass depends on the EoS, the nota¬ 
tion of a “relatively high” or “relatively low” binary mass 
is EoS-dependent as well. As a consequence, for Mtot = 
2.7 Mq all three types of mergers can occur depending 
on the EoS (see Fig. 5 in [22]). For EoSs which lead to 
compact NSs, the initial stars merge with a higher im¬ 
pact velocity (see Fig. 3 in CS), which is why the quasi- 
radial mode is strongly excited resulting in a pronounced 
/ 2-0 peak (Type I). For stiff EoSs, i.e. less compact stars, 
the merging proceeds more gently with a lower impact 
velocity, which suppresses the strong excitation of the 
quasi-radial mode. Such models favor the formation of 
pronounced spiral deformations, and thus the /spiral peak 
becomes particularly pronounced. The same reasoning ex¬ 
plains the occurence of the different merger types for a 
fixed EoS dependending on the total binary mass. 

We note that the classification of different merger sim¬ 
ulations in [22] is done based on the GW spectrum of the 
full signal (inspial and postmerger phase). Gonsidering the 
postmerger spectrum only (i.e. windowing the signal ap¬ 
propriately) the two secondary peaks are weaker, which 
suggests that they interfere with power from the inspiral 
signal and that their actual strength is weaker. While the 

® By windowing we mean the application of a window func¬ 
tion, e.g. a Tukey window, to the GW data, which can be used 
to exclude certain parts of the whole signal. 
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frequency of the inspiral signal may not even reach such 
high frequencies, some power may still be present at the 
frequencies of the secondary peaks because the inspiral sig¬ 
nal is finite. We also note that this classification scheme, 
the description of the different mechanisms and the follow¬ 
ing discussion of the frequency dependencies hold for sym¬ 
metric and mildy asymmetric binary mergers. The cases 
of strongly asymmetric mergers still need to be clarified 
as well as the nature of several other features which are 
even weaker than the secondary peaks. 

Finally, we mention that windowing the GW data has 
an impact in particular on the precise frequency of the 
/spiral peak. Choosing the postmerger spectrum only or 
even excluding the very early postmerger phase, the /spiral 
peak shifts to higher frequencies and its strength decreases. 
The latter is understandable because one excludes a part 
of the signal, when the mechanism generating the /spiral 
feature is operating. The frequency shift may not be unex¬ 
pected, given that the /spiral peak is produced by a highly 
dynamical feature. For instance, one may expect that the 
antipodal bulges orbit at slightly smaller radii and thus 
lead to slightly higher orbital frequencies. 

5.3 Frequency dependence of secondary peaks (and 
EoS constraints) 

The dependence of the dominant oscillation frequency /peak 
on the EoS was extensively discussed in Sect. |31 and es¬ 
pecially the potential for an accurate determination of NS 
radii (or, equivalently, NS compactnesses at fixed masses) 
was pointed out. The frequencies of the secondary peaks 
/spiral and / 2-0 follow a behavior very similar to that of 
/peak, i.e. for a fixed total binary mass the frequencies 
are higher for EoSs which yield smaller NS radii. Relat¬ 
ing the frequencies with the radii of nonrotating NSs (as 
in Sect. ID one finds similarly tight correlations for /spiral 
and somewhat less tight for / 2 _o (see na, see also m 
for a plot relating the secondary frequencies directly to 
the dominant frequency /peak)- The frequencies follow the 

order / 2 _o < fs piral ^ /peak- 

In principle, secondary frequencies may be employed 
for EoS constraints [1171111811^ . however, a few issues are 
worth being mentioned. The secondary peaks are weaker 
than the main peak even if one takes into account the 
better sensitivity of GW detectors at lower frequencies. 
Hence, secondary frequencies will be more difficult to mea¬ 
sure. In addition, the secondary peaks are broader in com¬ 
parison to /peak and often do not stand out clearly from 
the background. This will further impede an accurate de¬ 
termination of the secondary frequencies. Also, a detection 
of a single secondary peak will require additional informa¬ 
tion, e.g. from a measurement of /peak, to safely associate 
the secondary peak with either /spiral or / 2 - 0 - In compar¬ 
ison, a detection of the dominant peak frequency is more 
likely because of its strength, and, as discussed in Sect. El 
a single detection of /peak is sufficient to yield an accurate 
determination of NS radii. 

Given the similarity to /peak we do not further discuss 
the dependencies of the secondary peak frequencies here. 


but we refer to mss] for more details. Regarding the 
implications for EoS constraints we, however, note that 
our simulations in [22] do not confirm the existence of a 
mass-independent universal relation for the strongest sec¬ 
ondary peak as claimed in Instead, we find 

tight relations for the individual secondary frequencies for 
fixed total binary masses only. Note that there is no con¬ 
flict between the simulation data of mm and [22] . 
but different conclusions are explainable by the choice of 
the investigated binary setups. In fact, comparisons of the 
secondary and dominant frequencies for individual mod¬ 
els yield a good quantitative agreement. The relation pro¬ 
posed in [117U118] is built on simulations with a set of 6 
EoSs, but with different binary-mass ranges for each EoS. 
The choice of range of the binary masses, however, affects 
the distribution of GW frequencies. This is the reason why 
an EoS-dependent choice of binary mass ranges introduces 
a bias (depending on whether the mass range for a given 
EoS is relatively high or low compared to the average), and 
why a universal relation does not exist if the same range of 
binary masses is chosen for all EoSs as in . Apart from 
this, for a robust universality and mass-independence the 
mass range should comprise approximately the range that 
is expected from observed binaries, and not only a small 
variation of 0.2 Mq in Mtot- Finally, we note that a set 
of only 6 EoSs may not be sufficient to allow for robust 
conclusions about the spread in empirical relations and 
thus the quality of certain relations and their usability for 
EoS constraints, because models which may possibly lead 
to outliers are not included. Since a mass-independent re¬ 
lation of the secondary frequencies does not exist, the re¬ 
lation proposed in [11711118] cannot be employed for EoS 
constraints as proposed. However, as detailed in Sect. EJ 
detecting the weaker secondary GW peaks is not essen¬ 
tial for NS radius measurements if the stronger, dominant 
postmerger oscillation was measured. 


The understanding of the most prominent GW emission 
mechanisms as detailed in Sect. E| motivate us to set up 
an analytic model for the postmerger GW signal. This 
model may form the basis for GW templates to be used in 
matched filtering GW searches (see [35] for the potential 
of template-based searches in comparison to morphology- 
independent burst search algorithms). We report the model 
here by specifying the x-y component of the reduced quadrupole 
moment of the GW source, which can easily be interpreted 
as the cross polarization of the GW signal along the po¬ 
lar direction. The other components of the quadrupole 
moment can be deduced with the same set of parame¬ 
ters, from which the complete GW signal can be derived 

In [iTTirnsi no distinction between /spiral and / 2-0 was 
made, but their secondary frequency fi should correspond to 
the strongest secondary peak, which for most cases may be the 
/spiral feature. 


6 Analytic model for postmerger GW 
emission 
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(i.e. both polarizations in all emission directions). For in¬ 
stance, the plus polarization can be obtained by adding 
an appropriate phase shift. In |22| we identify three main 
mechanisms that produce the dominant postmerger GW 
emission. As also detailed in Sect.[5]this includes the domi¬ 
nant oscillation by the fundamental quadrupole mode, the 
coupling of this mode with the quasi-radial mode of the 
remnant and the transient emission from a spiral deforma¬ 
tion. All three mechanisms can be modelled by individual 
sine functions with given initial amplitude, initial phases 
and exponential damping behavior. Hence, our analytic 
model reads 

hx OC Qxy —Apeai^ exp ( (t to)/Tpeak) 
sin (27r/peak(l ^o) F ^peak) 

F Agpii-al exp ( (t to)/'^spiral) 

sin (27r/spii.al(^ ^o) F (ii'spiral) 

-f A 2-0 exp (-(< - to)/T 2 -o) 
sin (27r/2_o(t - to) F ^ 2 - 0 ), 


for t > to with a starting time to- The particular advan¬ 
tage of our model (for instance in comparison to the fits to 
GW signals described by [69]) is that the parameters are 
motivated by the underlying physical mechanisms. Hence, 
it is not difficult to choose appropriate values for these pa¬ 
rameters if one wants to reproduce a given GW signal, e.g. 
from a numerical simulation. The frequencies /peak, /spiral 
and / 2-0 can be chosen as found in the spectra. The am¬ 
plitudes and the damping time scales can be estimated 
from the time evolution of the GW signal and adjusted 
such that the correspoding peaks in the GW spectrum co¬ 
incide with numerically obtained spectra. Generally, the 
amplitude and the damping timescale of the fundamen¬ 
tal mode (“peak”) should be the largest. However, for the 
model discussed above (Fig. [5]), Agpirai may have a com¬ 
parable strength but a shorter damping time scale than 
the contribution from the fundamental mode. As argued 
in [^ , the antipodal bulges of the spiral deformation dis¬ 
appear after a few orbits, which suggests a damping time 
scale Tgpirai of the order of a few milliseconds. 

The physical background of the different parameters is 
also advantageous because the individual parameters can 
be constrained to certain ranges. For instance, the damp¬ 
ing time scales T 2 _o and Tspirai are always smaller than 
Aieak- The restricted range of certain parameters is an 
important property of our model because it signficantly 
reduces the parameter space and thus the computational 
costs in future applications as templates in matched fil¬ 
tering searches. The range and number of parameters (for 
instance the amplitudes and the frequencies) may be con¬ 
strained if information about the total binary mass (or at 
least the chirp mass) is available from the GW inspiral 
signal. In this context, we point out that for a fixed total 
binary mass and a fixed mass ratio the three frequencies 
are highly correlated mm- For symmetric binaries with 
Mtot = 2.7 Mq, /spiral and / 2-0 are tight functions of 
/peak, which can be well approximated by 


( 8 ) 



Fig. 12. GW signal of a 1.35-1.35 Mq merger with the DD2 
EoS (black). The analytic model GW signal is shown as dashed 
blue line (hardly visible since overlaid by the black curve) and 
is virtually indistinguishable from the actual signal. 


and 

/2-o = 1.0024/peak- 1.0798, (9) 

with the frequencies in kHz. Since these relations are accu¬ 
rate with deviations of typically only a few 10 Hz, /spiral 
and / 2-0 can be essentially eliminated from our model, 
which reduces the dimensionality of the parameter space. 
If this precision is not sufficient for template searches one 
can at least define a very narrow range of frequencies 
around the estimates of Eqs. ([8|) and Further corre¬ 
lations and constraints of the different parameters will be 
explored in future work. For an application of this model 
in a template-based GW search, which we leave to fu¬ 
ture work, possible refinements of the damping behavior 
or small frequency drifts may be considered as well. Also 
a certain rise time of the signal may be implemented or 
the model may be directly connected to an inspiral signal 
(see EHI). For strongly aymmetric binaries modifications 
may be necessary. 

We note that our model naturally leads to a time 
varying instantaneous GW frequency if the signal is in¬ 
terpreted as being produced by a single instantaneous fre¬ 
quency (see |5^lll8l[7i] l. The interpretation of a single in¬ 
stantaneous frequency, however, cannot be supported by 
the underlying physical mechanisms producing the GW 
signal (see Sect. [S] and Instead, there is evidence 

that several mechanisms with different frequencies con¬ 
tribute simultaneously to the signal. (See also the time- 
frequency map in m, which clearly shows distinct fre¬ 
quencies being simultaneously present and relatively sta¬ 
ble in time.) This feature also distinguishes our analytic 
model from the model in m , which describes only a sin¬ 
gle instantaneous orbital frequency which contributes to 
the GW signal and thus neglects the presence of the orbit¬ 
ing antipodal bulges generating the distinct /spiral peak. 

To exemplify the potential of our analytic model we 
show in Fig. [T2] the postmerger GW signal of a simula- 


/spiral = 0.8058/peak- 0.1895, 
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Fig. 13. Initial phase of the postmerger GW signal of a 1.35- 
1.35 Mq merger with the DD2 EoS (black). The analytic model 
GW signal is shown as dashed blue line. 



Fig. 14. Postmerger GW spectrum of a 1.35-1.35 Mq merger 
with the DD2 EoS (black) given by he« = ^(/) ' / with the 
Fourier transform of the waveform ■ The spectrum obtained 
from the analytic model in shown as dashed blue line. Thin 
dashed lines show the anticipated unity SNR sensitivity curves 
of Advanced LIGO [ST] (red) and of the Einstein Telescope | 66 | 
(black). 


tion with the DD2 EoS (black) together with the signal 
of our analytic model (dashed blue curve, Eq. Q) com¬ 
puted for a chosen set of parameters. As described above, 
the frequencies are taken from the GW spectrum and ap¬ 
propriate values are chosen for the remaining parameters. 
We stress that for this example the parameters were not 
determined by a htting procedure but simply by a visual 


trial-and-error comparison of the numerical waveform and 
the model signal. The very good match between the wave¬ 
forms in Fig. [T^ is obtained without extensive fine-tuning 
of the parameters, but only by a crude, physically mo¬ 
tivated choice of the parameters. An even better match 
can be achieved by an appropriate fitting procedure or 
more elaborate tuning of the parameters of the model 
(see [Sll). This is shown in Fig. [T31 where we employ a 
publicly available algorithir@ to generate an even better 
fit to the data. For clarity we show only the initial post¬ 
merger phase but note that the late phase is equally well 
reproduced as in Fig. |T21 In Fig. [TH we show the cor- 
respodning spectrum of the numerical waveform and of 
the analytic model, which apparently reproduces well the 
prominent features. Note that it is simple to include an 
additional feature in our analytic model to account for 
the peak at about 3.8 kHz, which is the coupling of the 
quadrupolar mode with the quasi-radial mode appearing 
at a frequency / 2 - 1-0 = /peak + /o- Since most of the pa¬ 
rameters of this additional feature are already determined 
by the current model (Eq. d?])), the consideration of this 
/ 2 - 1-0 contribution would not significantly complicate the 
model. However, given the weakness of this peak and the 
lower sensitivity of GW detectors at higher frequencies 
we do not further investigate this option. Future work 
should investigate the performance of our model in fitting 
other waveforms, especially for asymmetric binary merg¬ 
ers, and the applicability of the model for template-based 
GW searches. 

Finally, we note that for the /spiral feature in the model 
shown in Figs. [T^ and im we have chosen an amplitude of 
Aspirai = 0.085 and a damping time scale of Tspirai = 3 ms. 
An amplitude of this magnitude corresponds to two point 
particles of about 0.2 Mq to 0.3 Mq orbiting at roughly 
the surface of the remnant. This further substantiates the 
finding in [22] that the /spiral peak in the GW spectrum is 
generated by antipodal bulges which form during merging 
and which orbit around the inner remnant for only a few 
milliseconds. 


7 Impact of rotation 

Most NS merger simulations have focussed on binaries 
with irrotational velocity profile, i.e. binaries with NSs 
that do not rotate intrinsically. This is justified by the 
conclusion that the viscosity is not sufficient to enforce 
co-rotation of the binary components by tidal interactions 
during the inspiral phase isniiiH]- In addition, the intrinsic 
rotation of NSs in binaries is expected to be relatively slow 
since the stars cannot be spun up by accretion because of 
a missing donor star. In fact, the fastest known pulsar in a 
NS binary system has a rotation period of only 22 ms (see 
e.g. list in [85]), which is slow compared to the orbital pe¬ 
riod of the binary prior to merging, which is of the order of 

® We use the Covariance Matrix Adap¬ 
tation Evolution Strategy downloaded from 
https://www.lri.fr/~hansen/cmaes_inmatlab.html, which 
is described in [60]. 
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f [kHz] 


Fig. 15. GW spectra of a 1.35-1.35 Mq merger with the DD2 
EoS with different initial intrinsic NS spins./leff = h{f) ■ f with 
the Fourier transform of the waveform hx- 


roughly 2 ms. An intrinsic rotation period of 22 ms is also 
slow in the sense that the stellar structure is practically 
unaffected by the rotation, which only becomes important 
for spin periods below ~ 5 ms. Finally, one should bear in 
mind that most NSs in binaries rotate even more slowly 
than 22 ms, and that in the particular case of a NS spin of 
22 ms the rotation will be further reduced by pulsar spin 
down due to magnetic dipole radiation until the binary 
components merge on a time scale of roughly 100 Myrs. 

Here we are mostly interested in the impact of intrin¬ 
sic NS rotation on the GW signal, in particular on the 
dominant peak frequency of the GW spectrum, which can 
be employed for accurate NS radius constraints. With re¬ 
gard to the discussion in Sect. E] the most crucial ques¬ 
tion is whether intrinsic NS rotation can substantially al¬ 
ter /peak for expected rotation rates compared to models 
without intrinsic rotation. (The relations like the ones in 
Fig. 121 could still be used to determine NS radii but with 
some loss of precision if the intrinsic rotation cannot be 
measured from the GW inspiral phase (but see e.g. [22] 
for a discussion of spin measurements).) Gomparing irro- 
tational, co-rotating and counter-rotating spins (with re¬ 
spect to the orbital motion) only very little influence on 
/peak was observed by [94] . More recent studies in full gen¬ 
eral relativity found a somewhat larger effect of the order 
of 100 Hz for intrinsic rotation rates roughly comparable 
to 22 ms [29] . A comparable influence was reported by [74] 
but for roughly three times faster rotation rate, suggest¬ 
ing a somewhat weaker effect. However, these results and 
the observed impact of rotation may not be representative 
because the chosen binary systems are close to the thresh¬ 
old to prompt BH formation, and most of these simulated 
NS merger remnants experience the delayed collapse to a 
BH a few milliseconds after merging. For such binary se¬ 
tups the postmerger object undergoes strong changes in 
the stellar structure, and therefore the GW signal may 
be more sensitively affected by small changes, e.g. by a 


small difference in the angular momentum of the remnant. 
Moreover, the calculations of [29] employed an ideal-gas 
EoS, which may not represent a proper description of NS 
matter at all densities. m considered unequal-mass bi¬ 
naries with fast intrinsic rotation ( exceeding the rotation 
rate of the fastest known pulsar in binary systems) and 
found a shift of /peak by only a few 10 Hz (in addition one 
setup with a misalignment between orbital and intrinsic 
rotation was investigated). 

In order to assess the impact of intrinsic rotation for 
a more representative setup we performed calculations of 
1.35-1.35 Mq binaries with the DD2 EoS and varied the 
initial intrinsic rotation rate using spin periods of 19 ms, 
7 ms, 4 ms, and 2.2 ms. These simulations were then 
compared to the results from a calculation with irrota- 
tional velocity profile, which represents the usual setup 
used for deriving the results in Sects. [3] to [S] Our simula¬ 
tions were performed with the code described in [231[251[T21 
I18j , which imposes the conformal flatness condition on the 
spatial metric for solving the Einstein equations [7311128] . 
For simulations of irrotational binaries this approximation 
showed a very good quantitative agreement with fully gen¬ 
eral relativistic calculations (comparisons were reported in 
e.g. [mmEZl?]). 

Figure [T2] demonstrates that the impact of initial NS 
rotation on the dominant postmerger GW frequency is 
practically negligible. Note that all tested models with 
NS spins rotate faster than the fastest known pulsar in a 
binary. For the model with the very fast rotation period of 
about 2.2 ms the peak frequency shifts by roughly 50 Hz 
relative to the model without initial NS spin, whereas 
the calculation with an initial NS spin period of 19 ms is 
practically identical to the computation without NS spin. 
One observes that models with faster initial spin lead to 
slightly lower peak frequencies and somewhat lower peak 
heights. This behavior is understandable from the fact 
that the remnants in such models have slightly more angu¬ 
lar momentum. Hence, they are less compact and oscillate 
at lower frequencies. Also, the higher angular momentum 
slightly damps the excitation of the f-mode oscillation be¬ 
cause of the centrifugal barrier, which is why the height 
of the main postmerger peak is slightly reduced. This can 
also be seen by considering the minimum of the central 
lapse function in the first compression phase (see Fig. [TU] 
for the model without intrinsic NS rotation), which reveals 
that mergers with initial spin are not as much compressed 
during the plunge and during the first compression phase 
as the corresponding model without spin. 

Finally, the insensitivity of the peak frequency with re¬ 
spect to variations of the initial NS spin is understandable 
by considering the angular momentum of the remnant. 
Gompared to the model without NS spin, the remnant’s 
angular momentum is only slightly higher in models with 
spin because most of the angular momentum originates 
from the orbital motion of the inspiralling NSs. The con¬ 
tribution from the NS spins is only at the level of a few per 
cent (only for the model with a spin period of 2.2 ms the 
additional angular momentum slightly exceeds 10 per cent 
of the orbital angular momentum). Hence, differences in 
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the remnant’s structure are marginal, which explains why 
the peak frequencies are hardly affected by initial rotation 
even if the NSs spin relatively fast in comparison to mea¬ 
sured spins in NS binaries. Therefore, the consideration of 
intrinsic NS spins can be safely neglected for the purpose 
of measuring the NS radii via the dominant postmerger 
oscillation frequency /peak (Sect. [3]). 

The different models displayed in Fig.[TH]show that the 
secondary peaks may be affected by very fast intrinsic NS 
rotation. While the model with a rotation period of 19 ms 
does not exhibit important differences compared to the 
calculation with an irrotational velocity profile, the bina¬ 
ries with faster spinning NSs lead to a stronger /spiral fea¬ 
ture. This makes sense considering our explanations about 
the origin of this peak (Sect. O. Clearly, initial NS rota¬ 
tion favors the formation of antipodal bulges and thus 
leads to higher /spiral peaks. Future work should also con¬ 
sider cases with different spin orientations and magnitudes 
of the individual binary components. 

8 Two-families scenario 

Recently, it has been proposed by m that two different 
families of compact stars, ordinary hadronic NSs and ab¬ 
solutely stable strange stars (SSs), might coexist under 
the assumption that the strange matter hypothesis is cor¬ 
rect [CTll29| (see e.g. [M1IM1I25] for similar ideas and [77] 
for arguments against absolutely stable strange matter). 
Specifically, the authors argued for a soft hadronic EoS 
and a stiff quark matter EoS. In this scenario, low-mass 
compact stars are NSs and compact stars with masses 
above a certain threshold are absolutely stable SSs, i.e. 
objects entirely composed of quark matter except for a 
dynamically unimportant crust of nucleonic matter with 
densities below the neutron drip density. A conversion of 
NSs to SSs may be triggered by different processes, e.g. 
mass accretion onto a NS. The mass-radius diagram in 
Fig- mi illustrates such a jump from the NS branch to 
the SS branch. During the conversion (e.g. [W1I551IM1I1001 
[i^48| l the baryon number of the compact star is con¬ 
served, and the change in the binding energy leads to a 
different gravitational mass. While clearly speculative, the 
scenario of two families of compact stars represents a pos¬ 
sible solution to the hyperon puzzle, i.e. the scenario is 
compatible with a hadronic EoS which is strongly soft¬ 
ened by the occurrence of hyperons, and a maximum mass 
of compact stars above ^ 2 Mq. Here we do not per¬ 
form actual simulations which treat the conversion pro¬ 
cess and the detailed postmerger evolution, but based on 
theoretical arguments we discuss the observational signa¬ 
ture of such a scenario without further commenting on the 
likelihood of the scenario or considering possible counter¬ 
arguments. For instance, it still remains to be explored to 
which extent the existence of SSs could lead to a cosmic 
flux of strangelets, which would prevent the existence of 
NSs because seed strangelets would convert a NS into a 
SS [Bl IHSlI^ . For SS EoSs which are compatible with the 
lower observational bound of the maximum mass of nonro¬ 
tating compact stars, we expect that significant amounts 



Fig. 16. Mass-radius relation of a hadronic EoS (black curve) 
and an absolutely stable strange quark matter EoS (dashed 
blue curve) within the two-family scenario |46| . M is the grav¬ 
itational mass and R the circumferential radius (EoSs provided 
by Giuseppe Pagliara). The dashed line indicates the conver¬ 
sion from a hadronic star to a quark star for a fixed rest mass. 


of strange quark matter become gravitationally unbound 
during a merger of SSs m- Whether there exists a cosmic 
flux of strangelets, however, depends also on the fragmen¬ 
tation and survival probability of the unbound strange 
quark matter (see e. g. [55] for a discussion of a possible 
reconversion of quark matter nuggets to ordinary nucle¬ 
onic matter). We focus on the GW emission of mergers of 
compact stars within the described scenario. 

We consider only one exemplary model of the two- 
families scenario, i.e. one pair of nucleonic and quark mat¬ 
ter EoSs (see Fig. [Tn|). (Specifically, for the hadronic EoS 
we use the relativistic mean field model of 113, which in¬ 
cludes hyperons and delta resonances and which is based 
on the model of m with density-dependent coupling 
constants. We consider the model where the ratios be¬ 
tween the coupling constants of the mesons to the A iso¬ 
bars and the coupling constants of the mesons to the nu¬ 
cleons are chosen to be Xa-A = 1-15, x^j a = ^,XpA = 1- 
The EoS of absolutely stable strange quark matter adopts 
the model of m with the parameter X fixed to 3.5. The 
EoS tables were provided by G. Pagliara.) Clearly, vari¬ 
ations to this model are possible, but the case discussed 
here is sufficient for the sake of outlining the special obser¬ 
vational signatures of such a scenario and the qualitative 
behavior. We expect no qualitative differences for other 
realizations of the two-family scenario. We consider only 
the GW emission of the postmerger phase and assume 
that several GW observations of compact-star mergers are 
made with different total binary masses (see also [3T] for 
an earlier discussion of the GW signature of SS mergers). 
We restrict the discussion to symmetric mergers. 

For the particular model discussed here the transi¬ 
tion of a NS to a SS occurs for NSs more massive than 
^ 1.45 Mq (gravitational mass), which corresponds to a 
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Fig. 17. Dominant postmerger GW frequency /peak as a func¬ 
tion of the total binary mass for symmetric mergers with a two- 
family scenario [46] . For low binary masses the merger remnant 
is composed of hadronuc matter (black curve), whereas higher 
binary masses lead to the formation of a strange matter rem¬ 
nant with a lower peak frequency (dashed blue curve). The 
vertical dashed line marks a lower limit on the binary mass 
which is expected to yield a remnant that is stable against 
gravitational collapse (see text). 

central energy density of Pconv = 1-5 x 10^® g/cm^ (see 
Fig. [TT)1) . This means that for total binary masses below 
2.9 Mq two NSs will merge, since the mass of each binary 
component is below 1.45 Mq. This should be considered 
as the generic case given that many observed NS binaries 
have a total mass below ^ 2.9 Mq. Prior to merging the 
maximum density in the initial stars starts to decrease be¬ 
cause of the deformations in the very last inspiral phase. 
During merging the densities then increase again (see e.g. 
Fig. 3 in ED). The maximum density in the merger rem¬ 
nant shows some oscillations and typically the maximum 
density in the remnant is larger than the one of the initial 
stars. However, in particular for low-mass mergers with 
Mtot below ^2.4 Mq, the maximum densities do not sig¬ 
nificantly exceed the one of the inspiralling stars because 
of the strong rotational (and thermal) effects in the rem¬ 
nant. Hence, for lower total binary masses the densities 
in the remnants remain below pconv and are thus not suf¬ 
ficient to trigger the conversion of the NS matter to ab¬ 
solutely stable quark matter (here we neglect a possible 
influence of thermal effects on the conversion threshold). 
For larger Mtot the densities in the remnanst may exceed 
the central density of a nonrotating NS with 1.45 Mq, and 
thus the conversion of NS matter to quark matter will set 
in. We note that a first density maximum is reached di¬ 
rectly in the very first compression phase during merging. 
During the subsequent evolution of the remnant the den¬ 
sities are typically lower than this first maximum. Only 
after many ms, the redistribution of angular momentum 
in the remnant may lead to an increase of the maximum 
density above the first maximum. Hence, the conversion to 
absolutely stable strange quark matter is triggered right 


after merging or it occurs only on a longer time scale, after 
the remnant oscillations and thus the GW emission have 
ceased. In summary, this results in the following picture. 
Binaries below a certain binary mass threshold Mtot.conv 
lead to a merger remnant made of NS matter. For to¬ 
tal binary masses above Mtot.conv the remnant encounters 
conditions that trigger the immediate conversion of NS 
matter to quark matter right during merging. Here we as¬ 
sume that the conversion takes place on a very short time 
scale that is shorter or at least comparable to the dy¬ 
namical time scale of the remnant [MlllHj. Hence, a quark 
matter remnant forms essentially immediately, which has 
the same total rest mass (but a somewhat reduced gravi¬ 
tational mass compared to the corresponding NS merger 
remnant). 

By performing a simulation for a 1.2-1.2 Mq binary 
with the hadronic EoS, we find that the maximum den¬ 
sity after merging is only somewhat lower than the tran¬ 
sition density Pconv, which triggers the conversion of a NS 
to a SS in a static configuration. From this we estimate 
that the binary mass which triggers a conversion is only 
slightly higher than 2.4 Mq, and thus we adopt in the 
following Mtot.conv = 2.5 Mq. We stress that thermal ef¬ 
fects may actually trigger the conversion (during merg¬ 
ing) already at lower densities. Hence, Mtot.conv may be 
somewhat lower, which we neglect for simplicity because it 
would require a detailed microphysical calculation of the 
temperature-dependent threshold. 

Based on these remarks we can now estimate the post¬ 
merger GW emission of these binary mergers without ac¬ 
tually performing simulations, which are technically non¬ 
trivial because the burning of NS matter to quark matter 
needs to be treated in an appropriate way including for 
instance an implementation of turbulent combustion. For 
the model discussed here we find i?i.6 = 10.47 km for the 
hadronic EoS and i?i.6 = 14.53 km for the SS EoSs by 
computing the TOV solutions (see Fig. |Tn|). Using these 
properties we can estimate the peak frequency of the post¬ 
merger phase by employing the fits discussed in Sect. [31 
which relate i?i,6 and /peak- From the fits we can com¬ 
pute the expected peak frequency for total binary masses 
of 2.4 Mq, 2.7 Mq and 3.0 Mq for both EoSs. With these 
data we can then interpolate linearly to predict the GW 
peak frequency as a function of Mtot, where we distinguish 
the cases Mtot < Mtot.conv and Mtot > Mtot.conv In the 
former case we employ the estimates from the hadronic 
EoS, while for the latter case we use the peak frequency 
derived from i?i,6 = 14.53 km, i.e. the quark matter EoS. 
In Fig. |T7| one can read off which GW peak frequency is 
expected for a given total binary mass, which is measured 
by the inspiral GW signal (see Sect. |2]). It is clear that 
the formation of a SS remnant leads to dramatic changes 
of the postmerger GW signal compared to the GW emis¬ 
sion from a NS remnant, which occurs for low total bi¬ 
nary masses. Figure [T7| demonstrates that two families of 
compact stars result in a jump to lower peak frequencies 
at Mtot.conv Such a scenario can be clearly distinguished 
from only one family of compact stars, because in this case 
one expects a continuous increase of /peak with Mtot (see 
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Fig 1 in ESI)- Observationally, this requires the detection 
of several merger events with different total binary masses 
to probe whether or not a jump to lower peak frequencies 
occurs with increasing binary mass. 

The vertical dashed line in Fig. [T7] displays the thresh¬ 
old mass for prompt BH formation for the hadronic EoS 
estimated via Eq. Thus, at least up to this binary 
mass we expect the remnant to be stable independent of 
how fast the onset of the conversion to strange matter can 
stabilize the remnant. 

We note that also SS EoSs roughly follow the /peak — 
i?i.6 relation, which justifies the use of the relations dis¬ 
cussed above (see [ISlEI)- Given the small spread in the 
fits to the empirical data of the /peak — Ri.& relations, the 
prediction of a frequency jump is a very robust feature 
of the two-family scenario. We remark that the transition 
from hadronic to quark matter leads to a change in the 
gravitational mass, which we do not take into account. 
Since these changes are small, they have only a secondary 
impact on the quantitative behavior in Fig. 1171 

The results in Fig. [17] represent the realization of only 
one possible model of the two-family scenario. However, 
we expect a qualitatively similar effect for other EoS mod¬ 
els within the two-family scenario. In this context it is 
worth mentioning that for the particular model discussed 
here (i.e. the pair of hadronic and quark EoSs) the total 
binary mass Mtot.conv for formation of a SS remnant is 
relatively low compared to the expected mass range of NS 
binaries. As argued in Sect. |2j binaries with total masses 
below 2.4 Mq should be considered unlikely. Hence, the 
probability to detect a binary with Mtot < -^tot.conv uiay 
not be high. However, we point out that there may be a 
second way to discern the two-family scenario from the or¬ 
dinary model of only one family of compact stars, which 
may succeed with only one detection of a binary in the 
most likely mass range of about 2.7 Mq. GW detectors 
may measure finite-size effects during the last cycles of 
the inspiral . Therefore, 

if the EoS information revealed from the inspiral phase is 
in stark conflict with the EoS constraint from the post¬ 
merger phase, this could be a strong indication for the 
two-family scenario, too. Note that for binaries in the most 
likely mass range, the inspiral probes only the hadronic 
EoS as long as Mtot < 2 x 1.45 Mq for the model dis¬ 
cussed here. Thus, the GW signal from the inspiral alone 
may not be sufficient to exclude the two-family scenario 
if no binary merger with Mtot > 2 x 1.45 Mq is mea¬ 
sured. Only the consideration of both merger phases may 
be conclusive in this respect. Einally, we note that also the 
possibility of SS-NS mergers for asymmetric binaries with 
individual masses above and below Mtot,conv/2 should be 
considered. We leave this for future work. 


® The SS EoSs considered in [IS] describe bare SSs. While 
the inclusion of a nuclear crust would not have an impact on 
the peak frequency, which is mostly determined by the high- 
density regime of the EoS, the nuclear crust would lead to 
somewhat larger stellar radii (a few hundred meters) and thus 
would fulfill the /peak —7?r.6 relation even better than indicated 
in Fig. 11 of [T8| . 


9 Conclusions 

We summarize the main ideas to infer EoS and NS prop¬ 
erties from the postmerger GW emission of NS binaries 
and in particular we review the new findings of this work. 

— The minimum requirement for accurate constraints of 
NS properties from the GW emission of NS merger 
remnants is the measurement of the total mass of the 
binary via the GW inspiral signal. The knowledge of 
the binary mass ratio is not critical, since the rem¬ 
nant’s dominant oscillation frequency depends only weakly 
on the initial binary mass ratio. We illustrate that in 
the context of NS mergers the total binary mass can 

be estimated very well from the chirp mass, which will 
be measured with very high accuracy. 

— For a fixed total mass the dominant oscillation fre¬ 
quency of the NS merger remnant scales tightly with 
the radii of nonrotating NSs of chosen fiducial mass 
for different EoSs. This in turn allows to determine 
the NS radius from a measurement of the dominant 
postmerger GW frequency, /peak- The error is given 
by the maximum deviation from the empirical relation 
between the GW frequency and the NS radius found 
among all EoSs. Choosing a fiducial NS mass some¬ 
what higher than the mass of one of the merging NSs 
minimizes the deviations in the empirical relation and 
allows NS radius measurements with an accuracy of 
about 100 to 200 meters. Consequently, already a sin¬ 
gle detection of only the dominant peak frequency is al¬ 
ready sufficient for a tight constraint on the NS radius 
and thus the high-density EoS through the frequency- 
radius relation. 

— We investigate asymmetric and symmetric binaries of 
the same chirp mass, which can be measured very accu¬ 
rately during the inspiral phase. The frequency-radius 
relation for the combined data from symmetric and 
asymmetric binaries is very tight (with maximum de¬ 
viations below 300 meters). Therefore, even under the 
pessimistic assumption that only the chirp mass was 
measured accurately and no information on the binary 
mass ratio is available, the frequency-radius relation 
can be employed for accurate NS radius measurements. 

— For a broad range of binary masses (in the representa¬ 
tive range between 2.4 Mq and 3.0 Mq) a single rela¬ 
tion between the dominant postmerger GW frequency 
and the NS radius can be constructed by rescaling 
the frequency with the total binary mass. Deviations 
from this relation are of the order of 500 meters. Thus, 
for practical purposes, relations for fixed total binary 
masses may be more useful because they are tighter 
and the total binary mass will be known. 

— From our representative simulations we conclude that 
intrinsic rotation of the initial NSs has no significant 
impact on the dominant oscillation frequency for initial 
rotation rates as expected in NS binaries. Thus, intrin¬ 
sic NS rotation can be safely neglected for deriving NS 
properties and EoS constraints from GW signals of the 
postmerger phase. 

— We present three different possibilities to constrain the 
maximum mass of nonrotating NSs from the collapse 
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behavior of NS mergers or from the measurement of 
the dominant GW frequency. One of the key observa¬ 
tions is that the binary threshold mass for prompt BH 
formation is well described by pure TOV properties, 
specifically by the maximum mass of nonrotating NSs 
and by the compactness of the maximum-mass TOV 
configuration. Thus, the threshold mass, which can be 
observationally determined, depends in a specific way 
on pure EoS properties. Also other properties of the 
maximum-mass configuration of nonrotating NSs like 
its radius or the maximum central density can be ob¬ 
tained from measuring postmerger GW emission, thus 
probing the very high density regime of NSs. 

— Our analysis provides further evidence that the dom¬ 
inant oscillation mode of the merger remnant is the 
fundamental quadrupolar fluid mode. In addition, our 
calculations show that the merger remnant forms a 
single, self-gravitating object right from the time of 
its formation. Although there are overdensities visible 
in the evolution of the density in the equatorial plane, 
which appear as a rotating double-core structure, these 
features should not be considered as independent dy¬ 
namical features, but as tracers of the nonaxisymme- 
tries of a single star. 

— Besides the dominant GW emission frequency (dis¬ 
cussed above), we identify two distinct processes which 
produce prominent secondary peaks in the GW spec¬ 
trum. One GW peak is prodcued by the orbital motion 
of an outer antipodal spiral pattern that forms during 
merging and persists for several milliseconds. Another 
GW peak is generated by the coupling of the funda¬ 
mental quadrupolar mode with the quasi-radial mode. 
Depending on the exact system parameters (EoS and 
total binary mass) one or the other secondary feature 
is more prominent or, for some cases, both GW peaks 
can have comparable strengths. The presence and the 
prominence of the different features follows a clear be¬ 
havior and can thus be embedded in a classification 
scheme of the postmerger dynamics and GW emission. 

” We present details of an analytic model for describ¬ 
ing the structure of the postmerger GW signal that 
is based on the physical mechanisms producing the 
GW emission. A given postmerger GW signal can be 
well described by our analytic model. The model pa¬ 
rameters are physically motivated and thus specific, 
physically motivated bounds on the parameters can be 
imposed, which is advantageous for a template-based 
GW search by reducing the computational effort. We 
point out the existence of certain parameter correla¬ 
tions which further reduce the dimensionality of the 
problem, and thus our description may form the basis 
for future matched-filtering GW data analysis. 

— We finally explore a more speculative scenario of two 
families of compact stars. Within this scenario massive 
compact stars are made of absolutely stable strange 
matter, whereas low-mass compact stars are ordinary 
NSs. Observationally this scenario can be revealed by 
a strong discontinuity in the dependence of the dom¬ 
inant postmerger GW frequency on the total binary 


mass. For the falsification or verification of this sce¬ 
nario the postmerger GW signal may be crucial, since 
the inspiral phase may only test the hadronic regime 
for the expected range of binary masses. 

Overall, these conclusions show that the postmerger 
evolution is a highly interesting phase for understanding 
NS and EoS properties. In particular, it seems that certain 
characteristics are only accessible through the postmerger 
phase. Clearly, apart from advances in the modelling of the 
different aspects of this object, further studies of the data 
analysis capabilities of the advanced GW detectors are 
required. The highly rewarding prospects sketched in this 
work strongly motivate the construction of even more sen¬ 
sitive GW detectors, such as upgrades of Advanced LIGO 
and Advanced Virgo or the Einstein Telescope, e.g. [5^ 
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